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Abstract 

In this work we examine the possible existence of new physics beyond the standard model 
which could modify the branching fractions of the leptonic (mainly tauonic) decays of bot- 
tomonium vector resonances below the BB threshold. The decay width is factorized as the 
product of two pieces: a) the probability of an intermediate pseudoscalar color-singlet bb state 
(coupling to the dominant Fock state of the Upsilon via a magnetic dipole transition) and a 
soft (undetected) photon; b) the annihilation width of the bb pair into two leptons, mediated 
by a non-standard CP-odd Higgs boson of mass about 10 GeV, introducing a quadratic de- 
pendence on the lepton mass in the partial width. The process would be unwittingly ascribed 
to the T leptonic channel thereby (slightly) breaking lepton universality. A possible mixing 
of the pseudoscalar Higgs and bottomonium resonances is also considered. Finally, several 
experimental signatures to check out the validity of the conjecture are discussed. 
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1 Introduction 



The Standard Model (SM) has become nowadays the necessary reference to confront exper- 
imental data with theory: any possible discrepancy between them is commonly denoted as 
New Physics (NP), actually implying the need for some new assumptions or extensions of 
the basic physical postulates. In the SM, matter and gauge fields follow different statistics, 
the former being fermions and the latter bosons. As is well-known, there are important 
reasons to believe that this is quite unsatisfactory. One of the major motivations to extend 
the SM is to resolve the hierarchy and fine-tuning questions between the electroweak scale 
and the Planck scale. Supersymmetry is a very nice solution in this regard, since diagrams 
with superpartners exactly cancel the quadratic divergences of the SM diagrams. In par- 
ticular, it requires that the boson sector of the SM including the Higgs structure should 
be enlarged with new scalar fields. Until recently, supersymmetry was thought as the only 
possibility to solve the hierarchy problem, partly because of the lack of known alternatives. 
However, a new formulation for the electroweak symmetry breaking (dubbed "little Higgs" 
theories [1]) has recently emerged where cancellation of divergences occur, conversely to 
supersymmetry, between particles with the same statistics. The (initially massless) Higgs 
fields can be seen in this framework as Goldstone bosons, adquiring a mass and becoming 
pseudo-Goldstone bosons via explicit symmetry breaking at the electroweak scale, but still 
protected by an approximate global symmetry which would keep them relatively light. 

This is not the whole story, however. Another approach to solve the hierarchy problem 
aimed to the old idea on Kaluza-Klein extra dimensions, either in an ADD scenario [2, 3] 
or in a Randall-Sundrum model [4]. In both cases, the scalar sector would be increased 
by the presence of neutral bosons like scalar gravitons and radions (the latter associated 
to the quantum oscillations of the interbrane separation), eventually leading to measurable 
deviations from the SM [5, 6]. Moreover, let us cite another example of this kind of ex- 
tensions of the SM: the axion model, originally introduced in [7, 8] as a consequence of 
the spontaneous breaking of the global U(1)pq axial symmetry. Nowadays, the axion is 
a pseudoscalar field appearing in a variety of theories with different meanings including 
superstring theory, yielding sometimes a massless particle and others a massive one: in 
astrophysics the axion represents a good candidate of the cold dark matter component of 
the Universe. 

Although there are well established mass bounds (e.g. from LEP searches [9]) for the 
standard Higgs boson and the Minimal Supersymmetric Standard Model (MSSM), the 
situation can be different in more general scenarios where the tight constraints on the 
parameters of the theory do not apply, leaving still room for light Higgs bosons compatible 
with present data although the point is currently controversial (see, for example, [10, 11, 
12, 13]). As a suggestive example, let us mention that a possible mixing between Higgs 
bosons from a two doublet structure (which will deserve special attention in this paper) 
and resonances would alter the exclusion limits put by LEP data on light Higgs bosons 
since the branching fraction (BF) into tau pairs decreases considerably in this case [10, 14]. 
Moreover, in the so-called next to MSSM (NMSSM) a gauge singlet superfield is added 
to the MSSM spectrum [15, 16]: new CP-even and CP-odd Higgs bosons enter the game. 
For some choices of the parameters of this model, one can obtain very light pseudoscalar 
Higgs states evading the LEP constraints and whose detection might require some dedicated 
efforts at the LHC [17]. 
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The search for axions or light Higgs bosons in the decays of heavy resonances has sev- 
eral attractive features: first, the couplings of the former to fermions are proportional to 
their masses and therefore enhanced with respect to lighter mesons. Second, theoretical 
calculations are more reliable, notably with the advent of non-relativistic quantum chro- 
modynamics (NRQCD) [18, 19]. Indeed, intensive searches for a light Higgs-like boson (to 
be generically denoted by 4>° in this paper) have been performed according to the so-called 
Wilczek mechanism [20] in the radiative decay of vector heavy quarkonia like the Upsilon 
resonance (i.e. T— > 70°). To date, none of all these searches has been successful, but have 
provided valuable constraints on the mass values of light Higgs bosons [21]. 

Nevertheless, in this paper we develop the key ideas already presented in [22, 23] on a 
possible signal of NP based on the "apparent" breaking of lepton universality in bottomo- 
nium decays. Stricto sensu, lepton universality implies that the electroweak couplings to 
gauge fields of all charged lepton species should be the same; according to our interpreta- 
tion, the (would-be) dependence on the leptonic mass of the leptonic branching fractions 
(Bee ; ^ = e, /i, r) of T resonances below the BB threshold, if experimentally confirmed by 
forthcoming measurements, might be viewed as a hint of the existence of a quite light Higgs 
- of mass about 10 GeV - deserving a closer look. 

1.1 Two-Higgs Doublet Models 

In its minimal version, the SM requires a complex scalar weak-isospin doublet to sponta- 
neously break the electroweak gauge symmetry. As already commented, theories that try 
to resolve the hierarchy and fine-tuning problems imply the extension of the Higgs sector. 
Loosely speaking, the simplest way (i.e. adding the fewest number of arbitrary parameters) 
corresponds to assume an extra Higgs doublet, i.e. the Two-Higgs Doublet Model (2HDM) 
[21]. The Higgs content of this theory is the following: a charged pair (H^, two neutral 
CP-even scalars (h°, H°) and a neutral CP-odd scalar (A ) often referred as a pseudoscalar. 
Let us also note that diverse extended frameworks beyond the SM can lead to an effective 
theory at low energies equivalent to the 2HDM. On the other hand, there exist models 
with higher representations for the Higgs sector (e.g. Higgs triplets or the above-mentioned 
NMSSM) leading to more complicated structures [24]. 

Any two-doublet Higgs model has to cope with the potential problem of enhancing the 
flavor- changing neutral currents (FCNC). Several solutions have been proposed to overcome 
this serious difficulty In the Type-I 2HDM only one of the Higgs doublets couple to quarks 
and leptons and, since the process which diagonalizes the mass matrix of quarks equally can 
diagonalize the Higgs coupling, there is no flavor- changing vertex for the Higgs bosons at 
the end. (Note that in such case the Higgs coupling to b quarks is not enhanced.) Another 
extreme possibility to avoid FCNC's is based on the assumption that one Higgs doublet 
does not couple to fermions at all whereas the other Higgs' couples to fermions in the same 
way as in the minimal Higgs model. On the other hand, the Type-II 2HDM allows one of 
the Higgs doublet couple to the up quarks and leptons while the other Higgs doublet can 
couple to down-type quarks and leptons. This is the kind of model on which we shall focus 
in the following, excluding MSSM 1 since current limits rule out a very light pseudoscalar 

lr The Higgs sector of the MSSM can be viewed as a particular realization of a constrained Type II 2HDM with 
less parameters free. However, in this paper we are not considering the 2HDM as a low-energy approximation of 
the MSSM, but in more general grounds. 
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Higgs boson [9] as advocated along this work. Nevertheless, other alternative scenarios as 
those mentioned in the Introduction can not be discarded. 

Among other new parameters of the 2HDM, one of special phenomenological significance 
in this work is the ratio of the vacuum expectation values (v± t 2 of the Higgs down- and up- 
doublets respectively) usually denoted as tan (3 = v 2 /vi, where v\ + v\ — v 2 with v ~ 246 
GeV fixed by the W mass. Indeed, tan/3 governs the Yukawa couplings between Higgs 
bosons and fermions, thereby potentially enhancing the rate of processes forbidden by the 
SM, but allowed thanks to new contributions heralding the existence of NP. 

The layout of the paper is the following: in section 2 we tentatively introduce the 
hypothesis of a light non-standard Higgs boson which could modify the leptonic decay rate 
of T resonances; in section 3 we firstly apply time ordered second-order perturbation theory 
for a two-step process: prior photon radiation from the Upsilon leading to a pseudoscalar 
intermediate state followed by its annihilation into a lepton pair mediated by a CP-odd 
Higgs. Alternatively, we consider in subsection 3.2 the factorization of the decay width 
assuming the existence of Fock states in hadrons containing (ultra)soft photons as low- 
energy degrees of freedom in analogy to gluons in NRQCD. In sections 4 and 5, we focus 
on a 2HDM(II) model and the effects of the postulated NP contribution on the leptonic 
branching fraction are analyzed in the light of current experimental data: we conclude 
from a statistical test that lepton universality can be rejected at a 10% level of significance. 
Possible mixing between a pseudoscalar Higgs and r]b resonances is also considered, and its 
consequences on the hyperfine splitting between vector and pseudoscalar states. We finally 
gather technical details in three appendices at the end of the paper. 



2 Searching for a light Higgs-like boson in T leptonic decays 

The starting point of our considerations is the well-known Van Royen-Weisskopf formula 
[25] including the color factor 2 for the leptonic width of a vector quarkonium state without 
neglecting leptonic masses, 

rS = ^Ql x K( Xi ) (1) 

where a ~ 1/137 is the electromagnetic fine structure constant; mx denotes the mass of 
the vector particle (a T(nS) resonance in this particular case) and Qt> is the charge of the 
relevant (bottom) quark (1/3 in units of e); -R n (0) stands for the non-relativistic radial wave 
function of the bb bound state at the origin; finally, the "kinematic" factor K reads 

K(x e ) = (1 + 2x e )(l - Axe) 1 / 2 (2) 

where X£ = m 2 /m\. Leptonic masses are usually neglected in Eq.(l) (by setting K equal to 
unity) except for the decay into t + t~ pairs. Let us note that K(xe) is a decreasing function 
of xf. the higher leptonic mass the smaller decay rate. Such a^-dependence is quite weak 
for bottomonium and, consequently, we will assume that lepton universality implies the 
constancy of the width (1) for all lepton species. 

2 As is well-known, gluon exchange in the short range part of the quark-antiquark potential makes significant 
corrections to Eq.(l) [26], but without relevant consequences in our later discussion as we focus on relative differences 
between leptonic decay modes. 
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However, in this work we are conjecturing the existence of a light Higgs-like particle 
whose mass would be close to the T mass and which could show up in the cascade decay: 



T 



0°h z + n 



£ — e, /j,,r 



(3) 



Actually, this process may be seen as a continuum radiative transition that in principle 
permits the coupling of the bottom quark- ant iquark pair to a particle of variable mass 
and J PC : ++ , CT + , 1 ++ , 2 ++ ... (always positive charge conjugation). In this investigation 
we will confine our attention to the two first possibilities: a scalar or a pseudoscalar bo- 
son. In fact, intermediate bound states and not the continuum will play a leading role in 
the process, as we shall see. In the language of perturbation theory, a magnetic dipole 
transition (Ml) would yield at leading order a pseudoscalar bb state from the initial-state 
vector resonance, subsequently annihilating into a dilepton. Alternatively, there should be 
a certain probability that a pseudoscalar color-singlet bb system could exist in the Fock 
decomposition of the physical Upsilon state, as the light degrees of freedom would carry 
the remaining quantum numbers. 

Throughout this work, we will focus on vector T states of the bottomonium family below 
open flavor 3 , and the complete process (3) actually would be 



where 7 S denotes a soft (= unobserved) photon and bb[n] stands for those intermediate 
states of different quantum numbers collectively denoted by n, either on the continuum or 
as bound states. Note that 0° is not a real particle in the channel (4), conversely to the 
Wilczek mechanism [20], but a virtual state mediating the annihilation of a bb intermediate 
state into the final-state lepton pair. Hence, if the 0° mass were quite close to the T mass, 
the Higgs propagator could enhance significantly the width of whole process. In fact, the 
analysis performed by OPAL [10] using LEP data assuming a mixing between Higgs bosons 
and bottomonium resonances [14] does not permit to exclude a light Higgs of mass around 
10 GeV using reasonable tan/5 values in the 2HDM(II). Moreover, one should not dismiss 
other possible scenarios, as pointed out in the Introduction. 

Since the radiated photon would escape detection in our guess 4 , the NP channel (4) 
would be unwittingly ascribed to the leptonic decay mode of the Upsilon resonance, intro- 
ducing in its Bu a quadratic leptonic mass dependence opposite to that of Eq.(2) due to 
the Higgs coupling to fermions This is a cornerstone in our conjecture but likely of practical 
significance only for the r + r~ decay mode, where missing energy is experimentally required 
as one of the selection criteria [2 7] : events with photons of order 100 MeV would be included 
in the sample of tauonic decays, ultimately contributing to the measured leptonic BF. On 
the contrary, the electronic and muonic BF's would be affected to a much lesser extent, 
both because of: i) the smaller leptonic mass; ii) the experimental constraint on the re- 
constructed dilepton invariant mass, which restricts severely the energy of possible "lost" 
photons 5 . 

3 The Y(3S) state is excluded in the present analysis since only experimental data for the muonic channel [28] 
are currently available; see http://pdg.lbl.gov for regular updates. 

Experimental measurements of Bu include soft radiated photons which, however, have to be taken into account 
for a consistent definition of the leptonic widths [28], as we claim for the NP contribution advocated in this work. 

5 The leptonic mass squared with a final-state photon is given by m u = m\(l — 2E 1 /mr). Hence _E 7 is much 
more limited by invariant mass reconstruction of either electrons or muons than for tau's where such constraint is 
not applicable. I especially acknowledge N. Horwitz and the CLEO collaboration for correspondence in this regard. 



T -> 7a 66[ n ](-> 0° -> ■ £ = e,n,T 



(4) 
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3 Intermediate bb pseudoscalar states 

Along this section, we examine the role played by intermediate states in the process (4) 
according to two different schemes: firstly, we use time-ordered perturbation theory (TOPT) 
to deal with the formation, via an electromagnetic transition from the initial-state T, of 
a virtual bb state and its subsequent annihilation into a lepton pair. As an alternative 
approach, we rely on the separation between long- and short-distance physics following the 
main lines of a non-relativistic effective theory (NRET) like NRQCD [19] - albeit replacing a 
gluon by a photon in the usual Fock decomposition of hadronic bound states - and NRQED 
[29, 30]. Different results for the final widths in each approach come up, however; they 
are discussed in subsection 3.3. Finally, let us note that, despite we generically refer to 
the Upsilon (T) in our study, actually we are focusing on the T(1S) state because of more 
precise data on its tauonic BF (B TT ) w.r.t. the T(2S), and not yet available for the T(3S). 

3.1 Time-ordered perturbative calculation 

Let us write the amplitude for the process (4) using TOPT at lowest order: 

Tt ^ a ~ ^ m r -E n -k + ie (5) 

n 

The sum extends over all possible bb intermediate states with proper quantum numbers 
and energy E n , and k is the energy of the unobserved photon. Bottomonium states in a 
1 S'o configuration should dominate the sum, as we are facing the radiative decay of a T 
resonance. Intermediary states of higher angular momentum will be suppressed, since the 
associated electromagnetic transitions would involve higher multipole moments; Ml transi- 
tions to S-wave states with different principal quantum numbers will neither be considered, 
as they involve orthogonal wave functions in the non-relativistic limit. 

From expression (5) we also see that intermediate states with energies closer to the T 
mass are enhanced with respect to those of higher virtuality. In that sense, the continuum 
bb contribution, starting with a pair of BB mesons (E BB - ~ 10.56 GeV), is well above 
the T(1S) and T(2S) masses, and the main contributions to the decay would come from 
intermediate pseudoscalar ^-wave bound states, i.e. the rjb resonances, differing from the 
T resonances in virtue of the hyperfine structure. 

Squaring the amplitude (5) and including phase space integrations, the width of the 
process reads 



dk p-y 



(|&7.|H|T) 



x 



m T — E v — k + ie 
J dE ee Pe+e -\(£ + r\}l\r ]b )\ 2 2iT5(my-k-Ee i ) + ... (6) 



where the dots stand for other intermediate states contributing to the sum in (5). For 
explicit expressions of the particle densities p 7 , pe+e-, we refer the reader to [31]. The last 
integral amounts to the width of a _+ resonance with mass m v * = mx — k decaying into a 
pair of leptons, which we will denote as T v *^u. Section 4 will be devoted to the calculation 
of this annihilation via the proposed Higgs exchange. 
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The matrix element squared for the Ml transition between spin-triplet and spin-singlet 
S- wave states of quarkonium can be written in terms of the T(nS) — > 7 S rjl(nS) width after 
performing the integration over allowed photon states 6 with energy k, 

2, / Pl | ( » 7 . | H | T ) | 2 = r ?i, 7 . „ = (|y V (7) 

where we emphasize the off-shellness of pseudoscalar resonance by writing In the last 
step, we used a non-relativistic approximation to estimate the width Ty^ aT) , [32]. We 
shall return later to this formula in greater detail. Plugging the above result into Eq.(6) 
one obtains 



Kin / (.),, \-' 1 f A ,, / 



T- 



dk ; ~2 Tr;*^« (8) 



\2m b / 2tt J \m r - E Vb - k + ie 



Under the soft photon hypothesis, we have restricted the integration over photon energies 
to an upper limit A <C rn^, dictated by both experimental photon energy resolution and 
event selection criteria as pointed out in section 2. In that case we can safely retain the 
first term in the expansion of the rj b energy, E Vb = a/ m rib + k 2 ~ m Vb + k 2 /2m Vb + ... in 
the denominator of (8): 

T^ lll= l ^(Sl\ 2 \Lf h ik t r r„,_„l (9) 

'• 3 \2mJ \2nJ |m T - m„ + iTJ2 - k\ 2 % ") 

The instability of the rjb intermediate state has been taken into account by the substitution 
m Vb m Vb — ir vb /2, valid in the narrow-width case. A rough estimate of this width can 



be obtained through the pQCD relation 



Qj s {m b ) 
a s (m c ) 
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With reasonable values of the quark masses, the running strong coupling and the measured 
r T/c = 16 ± 3 MeV [28] , we get F Vb ~ 4 MeV. 

Consider now the quantity between curly braces in Eq.(9): 

I [ * (mx - mm -V + ry4 r '^» ^ I(A - A£ '»» < 10 » 

The mass difference m^ — m Vb = AEh s is the hyperfine splitting between the T(nS)—rjb(nS) 
partners. Different approaches suggest that AE hs lies in the interval ~ 35 — 150 MeV for 
n = 1 (see [33] for a compilation of theoretical results). If A > AEhs + r % , i.e. the 
range of energies for the unobserved photon is higher that the mass difference AEhs, the 
integration region in Eq.(10) comprises the full 7] b resonance contribution; taking A ~ AEh s 
would reduce the result by a half, and values of A below AE hs — T Vb would make X almost 
vanishing. Finally, a very large value of A would make (10) quadratically divergent (and 
then A can be seen as an ultraviolet regulator). 

6 Notice that there is no infrared singularity in the case of a magnetic dipole radiation. 
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Now, keeping track of our discussion at the end of section 2, those events with photons 
of energies up to several hundred MeV would pass the selection criteria employed in exper- 
iments measuring the tauonic BF of T resonances. Therefore, one can safely take in this 
case (£ = r) the upper limit A in the integrals of Eqs. (8-10) of the order of few hundreds 
of MeV, well above the quantity AE hs + T Vb , but always A <C m Vb . 

We thus evaluate the integral in Eq.(10) and take the limit of small T Vb , yielding 

Vb 

The off-shell width T v *^u has been transformed in this limit to the on-shell width r % ^«, 
as m rj * = m~£ — AEhs = In section 4 we will obtain an expression for this partial width 
dominated by the postulated Higgs boson. 

Finally, under the aforementioned assumption that the photon energy cutoff A is large 
enough, the partial width of the whole decay reduces to the factorized formula 



p mi F Vb ^ u 16a/Q 6 \ 2 3 T Vb ^ u 



where is now the on-shell Ml transition between real T(n) and r)b(n) states (i.e. 

k = AEh s ). Dividing both sides of Eq.(ll) by the Upsilon total width Tx, one gets the 
final BF as the product of two BF's, namely £>x^ 7s « = Br^ lsVb x Br) b ->u- We have thus 
proved that our approximation for the full process matches the corresponding expression to a 
cascade decay taking place through a r]£ intermediate state above threshold. Remarkably, no 
dependence on the A parameter is left, provided that the condition AEhs + F Vb < A <C m Vb 
is satisfied. 

As argued before, other possible intermediate contributions (i.e. BB continuum) are 
numerically less relevant due to their increasing mass difference (mi - m B |j) in the de- 
nominator of Eq.(5); in fact, one can check that the interference between 77^ and continuum 
contributions in Eq.(6), assuming that the matrix elements ( rjb 7 S | H | T ), ( BB 7 S | H | T ) 
have the same /c-dependence, gives a correction which is well below 1% for both the T(IS') 
and T(2S*) resonances using AE hs ~ 100 MeV and A ~ 120 MeV as reference values. Pure 
continuum contributions would be even more suppressed. 



3.2 Long- and short-distance factorization according to NRET 

Instead of supposing a final-state photon radiated by the Upsilon via a Ml transition as in 
the precedent section, we will now consider the soft 7 S incorporated as a dynamical photon 
into a Fock state of the resonance, in analogy to soft gluons in the framework of NRQCD. 
Admittedly, the strong interaction rules the hadronic dynamics but electromagnetism still 
has small but observable consequences, e.g. isospin breaking effects [34]. 

In a naive quark model, quarkonium is treated as a nonrelativistic bound state of a 
quark- ant iquark pair in a static color field which sets up an instantaneous confining poten- 
tial. Although this picture has been remarkably successful in accounting for the properties 
and phenomenology of heavy quarkonia, it overlooks gluons whose wavelengths are larger 
that the bound state size: dynamical gluons permit a meaningful Fock decomposition (in 
the Coulomb gauge) of physical states beyond the leading non-relativistic description with 
important consequences both in decays and production of heavy quarkonium. 
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Figure 1: (a) [upper panel]: Electromagnetic annihilation of a T(1S) resonance into a charged 
lepton pair through a virtual photon; (b) [lower panel]: Hypothetical annihilation of a bb^S^] 
state (existing either as a Fock component of the T resonance, or as a consequence of a Ml 
radiative transition) into a charged lepton pair through a Higgs-like particle (denoted by (p°). The 
vertical dotted line represents the separation between long-distance physics (on the l.h.s.) and 
short-distance physics (on the r.h.s.) corresponding to an arbitrary scale A which can be taken of 
the order of one hundred MeV, i.e. the soft photon energy. 

In particular, the bb system in a T vector resonance can exist in a configuration other 
than the dominant J p = 1~ since the light degrees of freedom mainly formed by gluons and 
light quark- ant iquark pairs can carry the remaining quantum numbers, albeit with a smaller 
probability. Moreover, it is arguable that soft photons should be included among those low- 
energy hadronic modes too, therefore allowing the heavy quark-antiquark system to be, for 
example, in a color-singlet, spin-singlet configuration, i.e. a 

\ bb l ls o } + Is) state 7 . Hence, 

such dynamical photons would participate in some decay modes of heavy quarkonium, in 
analogy to dynamical gluons. 

According to NRQCD, higher Fock states containing soft gluons indeed can participate 
actively in the decays of heavy quarkonium, as for example, its annihilation into light 
hadrons [19, 35]. On the one hand, the probabilities of possible Fock components are 
encoded in long-distance matrix elements; the heavy quark-antiquark annihilation itself, 
being a short-distance process, would be perturbatively calculable. On the other hand, 
details about the complicated nonperturbative hadronization of gluons into final-state light 
hadrons can be avoided in inclusive channels by assuming the hadronization probability 
equal to one. 

Likewise, one can consider Fock states containing dynamical photons, i.e. pertaining to 
the hadron during a long-time scale as compared with the short-time process represented 
in our case by the bb annihilation into a dilepton, of order \ jm^. Thus, a dynamical (quasi- 
real) photon, stemming from a | £»6[ 1 *S'q 1 ' > ] + 7 S ) Fock state, could end up (without need 
of hadronization) as a final-state photon - though undetected as postulated in this work. 
Obviously, important differences stand between dynamical gluons and photons so that it 
will be instructive to review several aspects of NRQCD relevant for our later development. 

7 We borrow this spectroscopic notation from NRQCD. Sometimes, the bb\^S^\ state will be denoted as rfc. 
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NRQCD is a low-energy effective theory for the strong interaction removing the un- 
wanted degrees of freedom associated to the heavy quark mass. The following hierarchy 
between hadronic scales is usually assumed in heavy quarkonium: itlq ^> itlqv ^> rngv 2 ~ 
Aqod, where uiq and v denotes the mass and relative velocity of the heavy quark 8 respec- 
tively, and Aqcd is the strong interaction scale. The bound state dynamics is chiefly domi- 
nated by the exchange of Coulombic gluons with four-momentum (E ~ mqv 2 , p ~ mv); soft 
gluons have four-momenta of order (rriQV, itiqv) and ultrasoft gluons of order (ijiqv 2 , rriQV 2 ). 
Likely the above hierarchy makes sense for bottomonium since v ~ 10 _1 ; then the energy 
of ultrasoft degrees of freedom turns out to be of order 500 MeV or less. In fact, the 
separation between energy scales is decisive in phenomenological applications of NRQCD. 

Bodwin, Braaten and Lepage [19] showed in a rigorous way that inclusive annihilation 
decays of heavy quarkonium can be factorized according to soft and hard processes: a) Long- 
distance physics is encoded in matrix elements, providing the probability for finding the 
heavy quark and antiquark in a certain configuration within the meson which is suitable for 
annihilation in each particular case; b) The short-distance annihilation of the QQ pair with 
given quantum numbers (color, spin, angular momentum and total angular momentum) 
which could be perturbatively calculated. Therefore, the decay width is written as 

T (H) = E f QQN( A ) ( H \°nW\H) (12) 

n 

where A is an ultraviolet cutoff of the effective theory, separating the high- and low-energy 
scales. The short-distance coefficient Tqq^(A) can be calculated as a perturbation series 
in oi s (2rriQ). The long-distance parameter (H\0^\H) determines the probability for the 
quarkonium to be in the [n]-configuration of the Fock decomposition and can be interpreted 
as an overlap between the QQ[n] state and the final hadronic state, requiring either a 
nonperturbative calculation (e.g. on the lattice) or the extraction from experimental data. 
Similar arguments apply to heavy quarkonium inclusive production (see, for example, [36] 
and references therein). 

In close analogy to the above procedure, let us also introduce an energy parameter A to 
separate short- from long-distance physics in the process under study in this work as shown 
in Fig.l. (One might identify numerically this A parameter with the upper limit of the 
integrals (8-10) in the precedent section.) Thus, a factorization of the decay width of the T 
into two pieces is applied: a) the probability V r {j]l^ s ) of the existence within the Upsilon 
of a | 66[ 1 *S'q 1 ' ) ] +7 s ) Fock state; b) the annihilation width T^^u into a lepton pair via Higgs 
boson exchange. Therefore, r T ^ 7s « can be written in this approach as the product 

Ty^ sU = V T {r 1 l ls )xT ril ^ u (13) 

This equation follows the spirit of the factorization given in Eq.(12). One important differ- 
ence, however, is that the long-distance quantity V r (f]^'~f s ) can be calculated perturbatively 
in QED using a quark potential model (for the initial and final state wave functions are in- 
volved) because of the smallness of the electromagnetic coupling a, in contrast to NRQCD 
matrix elements. On the other hand, the short-distance parameter T v *^u in Eq.(13) which 
we may identify with a particular ^qq^ coefficient in formula (12), can be calculated with 
the aid of the Feynman rules of the model under consideration, as we shall later see. 

8 We employ the symbol v to denote the relative three- velocity instead of v, usual in NRQCD, to avoid confusion 
with the vacuum expectation value of the standard Higgs boson. 
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3.2.1 Estimate of the probability P T (^7 S ) 

In NRQCD the probabilities for different Fock configurations of the heavy quark-antiquark 
pair in heavy quarkonium can be estimated according to the number and order of the chro- 
moelectric and chromomagnetic transitions induced by the interaction effective Lagrangian, 
needed to reach such states from the lowest configuration or viceversa. Let us point out 
that in our particular case we are considering a mixed situation where electromagnetic 
transitions occur between bound states of quarks which are mainly governed by the strong 
interaction dynamics. Moreover, another caveat is in order: we will compute the transition 
rate between on-shell states, whereas the pseudoscalar rjl state in (13) should be somewhat 
off-shell. Nevertheless, we will assume this off-shellness small on account of the low energy 
of the radiated photon. 

As commented at the beginning of section 3, we are focusing on the T(1S) resonance, 
mainly because of much more precise experimental data on its tauonic BF as compared to 
the T(2S'), as displayed in Table 1, and still missing for the T(3S). In addition, the larger 
number of possible intermediate pseudoscalar bound states for the two latter resonances 
would make more complicated the theoretical analysis as compared with the T(IS') state, 
where only the lowest bb^S^] Fock configuration (i.e. a single r] b (lS) state) should con- 
tribute to the final annihilation into a dilepton. Therefore, a textbook expression [32] has 
been employed to calculate the width corresponding to a transition between S-wave states, 
i.e. a direct Ml-transition between the T(IS') and the rjb(lS) resonances. 

The probability of the Fock state \rjl^f s ) "inside" a T(IS') resonance is then estimated 
as the ratio, 

pMl -i a pp, 

V^vlls) - - ± ^ AEl x \ Mr .,f ~ 1CT 4 (14) 

using mj = My/ 2 ~ 5 GeV, and the hyperfine mass splitting between the T and t]b states 
AEf ls ~ 50 MeV as a reference value 9 . The matrix element Air-. Vb is defined as 

POO 

Mr: Vb = air u r {r) j (kr/2) u Vb (r) 
Jo 

where Uij(r) represents the reduced radial wave function of the initial and final resonance 
respectively, and jo is the spherical Bessel function. We have made in (14) the reasonable 
approximation: M.r-. Vb ~ 1. Let us remark, however, that this parameter involves the wave 
functions of the T and i] b resonances, actually constituting a nonperturbative matrix ele- 
ment appearing in the long-distance part of the factorized width, in analogy to conventional 
NRQCD. 

In the absence of current experimental data, it is worth noting that the order-of- 
magnitude ~ 10 -4 of (14) agrees with more elaborated calculations. For example, Lahde 
[37] obtains for the partial width of the process T(IS') — > r/^lS)^ the value 7.7 eV for 
AE hs = 59 MeV, which corresponds to a BF of 1.5 x 1CT 4 in accordance with (14). 

9 In the language of a low-energy effective theory, those low-energy photons would be properly designed as 
ultrasoft in accordance with the energy scale hierarchy for bottomonium. The incorporation of quite higher energy 
photons to the Fock decomposition of the resonance would appear more problematic as the typical lifetime of such 
Fock states becomes comparable to the short time-scale of the annihilation process. 
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Let us now compare the partial widths of the T(IS') decays into three gluons and two 
gluons plus a photon by means of the ratio [26] , 

r(T(15') -> igg) 4af, „ a a s\ m-2 



'- wm^f) **^ 1 -™^ ~ 10 (15) 

at the energy scale /i 2 = m T . Actually, r can be crudely viewed as the relative probability 
(or suppression factor) of a |£»6[ 1 5'q 1 ' ) ] +7 s ) Fock state w.r.t. the corresponding 1 66[ 1 5 , Q 8 ' t ] + g) 
Fock state in the T resonance. On the other hand, when the gluon energy is of order 
m^v 2 or less, the probability for the latter colored Fock state scales as v 4 according to the 
NRQCD scaling- velocity rules [38, 19]. Thus, numerically V (r]l^ s ) ~ r x v 4 ~ 10~ 4 , since 
v 2 ~ 1CT 1 for bottomonium, showing the same order-of-magnitude as Eq.(14). Hence, we 
will assume that Eq.(14) provides a reasonable estimate of the probability for the existence 
of a r]l state within the Upsilon 10 as a function of AE hs . 

Combining the "master" formula (13) and equation (14), the NRET approach leads to 

r T ri lsrib r _ 116a/Q b \ 2 3 

for the final decay width, to be compared with Eq.(ll) obtained from TOPT. 



3.3 Discussion 

From Eqs.(ll) and (16) it becomes apparent that the calculations of the final width r r ^ Js u 
based upon TOPT or upon factorization a la NRET lead to different results. Indeed, the 
total width F Vb of the rji, resonance appears in the denominator of formula (11) after the 
approximations; instead, the total width Tt of the T resonance appears in formula (16). 
Since one expects Tx <C T % 11 , the decay rate for the whole channel T — > ^stt stemming 
from the NRET factorization turns out to be much larger than the decay rate obtained 
from TOPT. In other words, both approaches are not dual each other. 

Actually, the fact that one gets different expressions and hence dissimilar numerical 
values for the width r T ^ 7s « in both methods is not contradictory in itself for they are 
based on distinct physical assumptions. The time-ordered perturbative scheme with r/l 
production above threshold essentially implies a cascade decay, i.e. photon emission and 
subsequent annihilation of the intermediate hadronic state. Only under this hypothesis and 
the narrow width approximation, the factorization given by Eq.(ll) is justified. On the 
other hand, the factorization a la NRET postulated in Eq.(16) assumes that a pseudoscalar 
color-singlet bb state exists inside the T as a Fock state. Both starting points are different 
and the results stemming from each framework need not coincide 12 . 

10 One might also consider in our conjecture a bfe[ 1 S' ( j 1 ' > ] + 3;;} Fock state yielding an unobserved hadronic system in 
the final state (2n,p, ...). However, this contribution would be very much suppressed because of the large virtuality 
of the intermediate rfe state since certainly AEhs should be quite smaller than the hadronic invariant mass. 

11 Because the hadronic decay of a pseudoscalar resonance via the annihilation of the QQ pair into two gluons 
should proceed at a higher rate than the corresponding decay channel of the vector resonance via three gluons. 

12 Let us make a pedagogical analogy with radioactive nuclides: the factorization given in expression (11) amounts 
to the product of branching fractions in a cascade decay; conversely formula (16) corresponds to the coexistence of 
a radioactive nuclide in some proportion with a stable isotope in nature. The decay rate of a sample of this element 
would be given by the fraction of the radioactive isotope (i.e. the probability to find a radioactive atom in the 
sample, analogous to P r (r/^js)) multiplied by its decay rate. 
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In fact, equivalent situations can be found, e.g. in inclusive hadropro duct ion of heavy 
quarkonium at large transverse momentum [39], when comparing the color-singlet frag- 
mentation mechanism (where all three perturbative final-state gluons are attached to the 
hard interaction Feynman diagram) versus the color-octet mechanism (where the emission 
of two soft gluons from a nonperturbative colored state takes place on a long-time scale). 
Another example where higher Fock states can compete with perturbative calculations is 
the explanation given by Braaten and Chen of the long-standing p — n puzzle of J / ip and 
ip f decays [40]. In general, their decays into light hadrons should occur by the annihilation 
of the cc pair into three gluons but the discrepancy with experiment is about two orders 
of magnitude in this channel. Instead, they argue that the process is dominated by the 
higher color-octet contribution ccpS^], thereby annihilating into a light quark pair via 
cc — > g* — > qq. The suppression of this decay mode for the ip f is attributed to a dynam- 
ical effect which cancels the cc wavefunction at the origin. (Alternatively, Brodsky and 
Karliner [41] suggested that the decay into pjr should proceed through the intrinsic charm 
components of the light mesons.) 

Hereafter, we will adopt the NRET factorization as given by the master equation (13) 
and the numerical estimates will be based on Eq.(16). The low-energy regime of the (ul- 
tra)soft photons provides confidence in this approach 13 . 



4 Effects of a light neutral Higgs on the leptonic decay width 

Following a general scheme, fermions are supposed to couple to the 0° Higgs field according 
to a Yukawa interaction term in the effective Lagrangian, 

A2 = ^W(*75)/ (17) 

where £y denotes a factor depending on the type of the Higgs boson and the specific theory 
under consideration, which could enhance the coupling with a fermion (quark or lepton) of 
type "/" and therefore plays a crucial role in our conjecture. In particular, 0° couples to the 
final-state leptons proportionally to their masses, ultimately required because of spin-flip 
in the interaction of a fermion with a (pseudo) scalar, thereby providing an experimental 
signature for checking the existence of a light Higgs in our study. Lastly, note that the 275 
matrix stands only in the case of a pseudoscalar 0° field. 

In this paper we are tentatively assuming that the mass of the light Higgs sought stands 
close to the T resonances below BB production: m^o < 2m b . As will be argued from current 
experimental data in the next section, we suppose specifically that m^o lies somewhere 
between the T(lS) and T(2S) masses, i.e. 

ra-r(is) < "ty> < m r(2 s) (18) 

Now, we define the mass difference: Am = jm^o — m Vh \, where rj b denotes either a IS* or 
a 2S state. Accepting for simplicity that the Higgs boson stands halfway between the mass 
values of both resonances, we set Am ~ 0.25 GeV for an order-of-magnitude calculation. 

13 On the other hand, in using Eq.(ll) instead of Eq.(16) large values of tan/3 are required at the end of the 
calculation to account for the leptonic BF rise with the lepton mass, leading to a ijb width exceendingly large 
in contradiction with the narrow width approximation made along the way to get (11). On the contrary, no 
inconsistencies of that kind arise when using the NRET approach to interpret the conjecture made in this work. 
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Hence we write approximately for the scalar tree-level propagator of the 0° particle in 
the process (4) entering in the evaluation of T^^u, 

1 1 1 
~ ~ (19) 

« - m io) 2 + s° r J° ( m i- m i«y 4 <( Am ) 2 

where the total width of the Higgs boson T^o has been neglected assuming that (Am) 2 3> 
r| - We will make a consistency check of this point in subsection 5.1 using numerical values. 

Performing a comparison between the widths of both leptonic decay processes (i.e. T Vb ^u 
versus r^™^), one concludes with the aid of Eqs. (17-19) and (1) that 

^ 3mfmf(l - Ax^\RM\\lQ ^ 3gg mjm\ 1 (em) 

27r 2 (m2 6 -mj ) 2 v 4 ~ 32vr 2 Q 2 a 2 t; 4 Am 2 i; 4 1 + 2x £ *c—*u u ) 

Above we used the non-relativistic approximation (more precisely the static approxima- 
tion) when assuming null relative momentum of heavy quarks inside quarkonium, and the 

same wave function at the origin for both the T vector state and the ftftpSg 1 ^] intermediate 
bound state on account of heavy-quark spin symmetry [19]. Then, the decay amplitude 
squared of the pseudoscalar state into a J PC = ++ (CP-even) Higgs vanishes and only a 
CT + (CP-odd) Higgs would couple to pseudoscalar quarkonium in this limit. Therefore, the 
Higgs boson hunted in this way should be properly denoted by A and, consequently, this 
notation instead of the generic 0° will be employed for it from now on. 

4.1 Modification of the leptonic BF due to a light CP-odd Higgs contribution 

The BF for channel (4) can be readily obtained inserting Eq.(20) into Eq.(16) and afterwards 
dividing by Tx, as 



>7s« — 



8vr 2 a(l + 2x £ )T r Am 2 v 4 _ 



x B u (21) 



so one can compare the relative rates by means of the following ratio 



~ B„ 



£Ue m b AE hs 



2 



x 1 (22) 
_87r 2 a(l + 2xi)T^v A \ Am 2 

where we are assuming in the denominator that the main contribution to the leptonic chan- 
nel comes from the photon-exchange graph of Fig. 1(a). Let us point out once again that, 
since the 7 S remains undetected, the NP contribution would be experimentally ascribed to 
the leptonic channel of the T resonance. Thus, the ratio (22) represents the fraction of 
leptonic decays mediated by a CP-odd Higgs, ultimately responsible for the breaking of 
leptonic universality due to the quadratic mass term m 2 . 

Now, to facilitate the comparison of our results with other searches for Higgs bosons, we 
identify in the following the factor with the 2HDM (Type II) parameter for the universal 
down-type fermion coupling to a CP-odd Higgs, i.e. £b = & = tan/3 [21]. Inserting 
numerical values into (22) and keeping the leading term in m 2 , one gets the interval 

K ~ (1.5-10 -7 - 1.2-10" 5 ) x tan 4 /3 x m\ (23) 

where use was made of the approximation m^o ~ 2m& ~ 10 GeV, and the broad range 
35 — 150 MeV for the possible hyperfme mass diference AEh s [33]; me is expressed in GeV. 
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Table 1: Measured leptonic BF's (B u ) and error bars (a e ) in %, of T(15) and T(2S) (from [28]). 



channel: 


e + e 






T(1S) 


2.38 ±0.11 


2.48 ±0.06 


2.67 ±0.16 


T(2S) 


1.34 ±0.20 


1.31 ±0.21 


1.7 ± 1.6 



4.2 Possible A — i] b mixing 

Long time ago, the authors of references [42, 43] pointed out the possibility of mixing 
between a light Higgs (either a CP-even or a CP-odd boson) and bottomonium resonances 
(scalar or pseudoscalar, respectively). Later, Drees and Hikasa [14] made an exhaustive 
analysis of the phenomenological consequences of the mixing on the properties of both 
resonances and Higgs bosons. In view of new and forthcoming data on the bottomonium 
sector from B factories [44, 45], we are particularly interested to apply those ideas looking 
for experimental signatures to provide an additional check on our conjecture of a light 
CP-odd Higgs particle. 

On the one hand, the mixing can enhance notably the gg decay mode of the Higgs 
boson [14], ultimately increasing its total decay width. A net effect would be an important 
decrease of the Higgs tauonic BF (when the A mass ranges from 9.4 GeV to 11.0 GeV, 
especially if it lies close to the t/& masses [14]). An exciting experimental consequence arises 
in the search for Higgs particles carried out at LEP: higher tan/5 values are allowed than 
those upper bounds derived from the analysis without considering the mixing [10, 11, 12]. 

Interestingly, the mixing with a CP-odd Higgs could also modify the properties of the 
pseudoscalar resonances, i.e. the r] b states. Thus, even for moderate tan/3 one might expect 
a disagreement between forthcoming experimental measurements of the hyperfine splitting 
mx ^ m r lb i and theoretical predictions based on potential models, lattice NRQCD or pQCD 
[33]. We will come back to this discussion in our numerical analysis of subsection 5.1. 

5 Lepton universality breaking? 

Let us confront our predictions based on the existence of a CP-odd Higgs boson with 
experimental results on T leptonic decays [28] summarized in Table 1. Indeed, current data 
show a slight rise of the decay rate with the lepton mass when comparing the r + r~ decay 
mode with the other two (e + e~ and /x + /i~) modes. However, error bars (cr^) are too large 
(especially in the T(2S) case) to permit a thorough check of the lepton mass dependence 
as expressed in (23). Nevertheless, we have applied a hypothesis test (see appendix A) to 
check lepton universality using the T(IS') and T(2S) data displayed in Table 1. 

The null hypothesis (i.e. lepton universality) is compared against the alternative hy- 
pothesis stemming from the Higgs contribution predicting a larger (and positive) value of 
the measured mean of the £>«'s differences. Thereby, we conclude that lepton universality 
can be rejected at a 10% level of significance. As a cornerstone of this work, such slight 
but measurable variation of the leptonic decay rate (by a (9(10)% factor) from the elec- 
tronic/muonic channel to the tauonic channel can be interpreted theoretically according to 
the 2HDM upon a reasonable choice of its parameters (e.g. tan (3) as we shall see below. 
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5.1 Values of tan/5, A — r\ h mixing and discussion 

In order to explain the rise of the tauonic BF by a ~ 10% factor w.r.t. the electronic/muonic 
decay modes, one obtains from Eq.(23) that tan/5 should roughly lie over the range: 

7<tan/3<21 (24) 

depending on the value of AEh s , namely from 150 MeV down to 35 MeV, whose limits 
remain somewhat arbitrary however. The partial width for the tauonic decay mode of the 
T(1S) mediated by the CP-odd Higgs turns out to be r T ^ 7srr ~ 140 eV. 

A caveat is in order: the above interval is purely indicative since there are several 
sources of uncertainty in its calculation, like the actual mass of the hypothetical Higgs and 
not merely the guess made in Eq.(18) or the crude estimate of the probability V r (i] b *'-f s ). 
In fact, higher values of tan/3 were obtained in our earlier work [22] because lower photon 
energies were used (i.e. between 10 and 50 MeV). Actually, letting V r (rj b *'y s ) vary, the range 
given in Eq.(24) changes accordingly and somewhat higher values cannot be ruled out at 
all. In sum, our calculations are only approximate and we cannot claim a well-defined 
interval for tan f3 but just an indication on the values needed to interpret a possible lepton 
universality breakdown according to our hypothesis 14 . Nevertheless, we perform below a 
consistency check of (24) concerning several partial widths of the rj b and A particles. 

Firstly, let us insert the values of tan/3 given by Eq.(24) into Eq.(20) to compute r % _ >rr ; 
notice that a high value of tan (3 might yield a large partial width for the r)b resonance, as 
compared with the expectation T Vb ~ 4 MeV obtained in section 3.1. In fact, using the 
interval given in (24) one gets T Vb ^ TT varying from 44 keV up to 3.56 MeV. Therefore, taking 
into account the NP contribution to the total decay rate, the Higgs-mediated tauonic BF 
of the rjb resonance should stay over the range 0.01 < r Vb ^ TT /T rib < 0.5. 

On the other hand, the decay width of a CP-odd Higgs boson into a tauonic or a cc pair 
in the 2HDM(II) can be obtained, respectively, from the expressions: [14] 

r(A°^T+T-) ^ ^^,.(1-^ (25) 

T(A°^cc) ~ ^l 00 ^ m A0 (l-4x c )^ (26) 

where x T = m 2 ./m A0 and x c = m 2 /m 2 A0 . Below open bottom production, even for moderate 
tan (3, the A decay mode would be dominated by the tauonic channel, i.e. 

r A o ~ T(A° -> t + t-) — 1-10 MeV (27) 

Thus we can confirm the validity of the approximation made in Eq.(19) for the Higgs 
propagator, i.e. Am 2 ^> T A0 (where we tentatively set Am = 250 MeV). 

As commented in subsection 4.2, there is another interesting consequence of our conjec- 
ture related to bottomonium spectroscopy due to the mixing between a CP-odd Higgs and 
rjb states. (In appendix B we introduce the notation and basic formulae.) Indeed, using 

14 It is also worthwhile to remark that the range in (24) is compatible with the lowest values of tan f3 needed to 
interpret the g — 2 muon anomaly in terms of a light CP-odd Higgs resulting from a two-loop calculation [11, 12, 46]. 
At present, there is a discrepancy (3.0a) between the theoretical value and the experimental result based on e + e~ 
data, but only 1.0<t when r data are used [47]. Hence the situation is still unclear to claim for new physics beyond 
the SM from the g — 2 analysis alone. 
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the values of tan/5 from Eq.(24) the mixing parameter defined in Eq.(B.2) turns out to be 
5m 2 ~ 1.0 — 3.1 GeV 2 . Therefore, such r^ — A mixing could induce an observable mass shift 
of the physical rjb states which would eventually increase the hyperfine splitting between 
pseudoscalar and vector resonances w.r.t. a variety of calculations within the SM [33]. 

Let us now write the masses of the mixed (physical) states as a function of the masses 
of the unmixed states (marked by a subscript '0', i.e. r] b0 and Aq) with the aid of the 
expression derived from Eq.(B.5) for narrow states and m Vb0 ~ m A o, 

K,a° - 1«°+<J ± ^[Ko-< ) 2 + 4(W) 2 ] 1/2 (28) 

Taking as a particular case tan/3 = 16 (5m 2 ~ 2.3 GeV 2 ), m Vb0 = 9.4 GeV and m A o = 9.5 
GeV, we get approximately m Vb ~ 9.32 GeV and m A o ~ 9.58 GeV, compatible with our 
tentative hypothesis on the Higgs mass (i.e. Am = \m A o — m Vb \ ~ 0.25 GeV) and the 
experimental mass 15 so far measured for the rjb meson [33, 28]. 

6 Summary 

In this paper we have interpreted a possible breakdown of lepton universality in T leptonic 
decays (suggested by current experimental data at a 10% level of significance) in terms of a 
non-standard CP-odd Higgs boson of mass around 10 GeV, thereby introducing a quadratic 
dependence on the leptonic mass in the corresponding BF's. Higher-order corrections within 
the SM (involving the one-loop r/b decay into two photons, as estimated in appendix C) fail 
by far to explain this effect. 

The existence of a CP-odd Higgs of mass about 10 GeV mixing with pseudoscalar bb 
resonances should display very clean experimental signatures and therefore could be easily 
tested with present facilities: 

• "Apparent" breaking of lepton universality when comparing the BF of the t + t~ decay 
mode on the one hand, and the BF's of the electronic and muonic modes on the other. 
Experimental hints of this possible signature triggered this work [22]. 

• Presence of monoenergetic photons with energy of order 100 MeV (hence above detec- 
tion threshold) in those events mediated by the CP-odd Higgs boson (estimated about 
10% of all T tauonic decays). This observation would eventually become a convincing 
evidence of our working hypothesis. 

• The T — r/b hyperfine splitting larger than SM expectations, caused by the A — rjb 
mixing. Also a rather large total width of the r/ b resonance due to the NP channel 
(especially for the higher values of tan/3 shown in (24)). 

Although we have focused on a light CP-odd Higgs boson according to a 2HDM(II) for nu- 
merical computations, the main conclusions may be extended to other pseudoscalar Higgs- 
like particles with analogous phenomenological features, as outlined in the Introduction. We 
thus stress the relevance for checking our conjecture of new measurements of spectroscopy 
and leptonic decays of the Upsilon family below BB threshold in B factories (BaBar [44], 
Belle [45]) and CLEO [48]. 

15 The measured mass [33] of the observed r/b (IS 1 ) state is 9300±20±20 MeV, indeed slightly smaller than different 
calculations of the hyperfine splitting. This measurement based on a single event needs confirmation however [48]. 
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Appendices 

A Lepton universality breaking: hypothesis testing 

In Table 2 we show the differences between the branching fractions of the leptonic channels 
defined as A u > = Byy — B u £ = e, fi, r) for both T(IS') and T(2S) resonances, obtained 
from Table 1 in the main text. Dividing them by their respective experimental errors 
o~w — a/o"| + o~J, (see also Table 1) one gets the ratios A^/aw. Only four of these quantities 
can be considered as independent. Moreover, in view of the small difference between the 
electron and the muon masses as compared with the tau mass, we will base our analysis 
on the comparison between the electron and the muon decay modes on one side, versus the 
tauonic mode on the other side 16 . 

In this analysis, we are especially interested in the alternative hypothesis based on the 
existence of a light Higgs boson enhancing the decay rate as a growing function of the 
leptonic squared mass, in opposition to the kinematic factor (2). Therefore, the region of 
rejection for our statistical test should lie only on one side (or tail) of the A^j / aw variable 
distribution (i.e. positive values if my > me), in particular above a preassigned critical 
value [49, 50]. In other words, we have performed a one-tailed test [49, 50] using the sample 
consisting of the four independent BF differences between the electronic and the muonic 
channels versus the tauonic decay mode as explained above (i.e. A eT /o~ eT , A MT /o" MT ). For 
the sake of simplicity, we will assume that such differences follow a normal probability 
distribution with a mean of 0.7775, obtained from the four independent values 

(£ = e,fi). 

Next, let us define the test statistic: T = (Ae T /ae T )xy/N ~ 1.555, where N = A stands 
for the number of independent points. Indeed note that we are dealing with a Gaussian of 
unity variance after dividing all differences by their respective errors. Now, we will choose 
the critical value to be ~ 1.288 corresponding to a significance level of 10% in the test. 
Lepton universality plays the role of the null hypothesis predicting a mean zero (or slightly 
less), against the alternative (composite) hypothesis stemming from the postulated Higgs 
contribution predicting a mean value larger than zero 17 . Since experimental data imply 
that T > 1.288, we can reject the lepton universality hypothesis at a 10% level of significance 
18 . Certainly, this result alone is not statistically significant enough to make any serious 

16 In a prior work [22] the muonic and tauonic modes were confronted with the electronic mode. Our final 
conclusion remains the same as before. 

17 We are facing a situation where the null hypothesis is simple while the alternative is composite but could be 
regarded as an aggregate of hypotheses [50]: we are assuming normal distributions, with unit variance and mean 
fio = for the null hypothesis, while fj, a > l^o = for the alternative complex hypothesis. A significance level of 10% 
means that the null hypothesis will be rejected if the measured mean value of the A ee i /a ee i differences is greater 
than w 1.288/\/iV, where N denotes the total number of points. This condition is equivalent to require that the 
test statistic T defined above should be greater than 1.288. 

18 Let us recall that the significance level (or error of the first kind) represents the percentage of all decisions such 
that the null hypothesis was rejected when it should, in fact, have been accepted [49, 50]. 
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claim about the rejection of the lepton universality hypothesis in this particular process, 
but points out the interest to investigate further the alternative hypothesis stemming from 
our conjecture on the existence of a light Higgs. 

Table 2: All six differences (obtained from Table 1 in the main text) between the leptonic BF's 
(in %) of T(15) and T(25) resonances separately. Subscript It' denotes the difference between 
channels into 11 and ft' lepton pairs respectively, i.e. A^/ = Bin> — the aw values were 

obtained from Table 1 by summing the error bars in quadrature, i.e. ou' = \J &f + cr 2 , . Only two 
A.££>/(T£gi values for each resonance can be considered as truly independent, amounting altogether 
to a total number of four independent points. 



channels 








r(is) eii 


0.1 


0.125 


+0.8 


t(isv 


0.19 


0.17 


+1.12 


r(is) eT 


0.29 


0.19 


+1.53 




-0.03 


0.29 


-0.10 


T(2S) I1T 


0.39 


1.61 


+0.24 


T(2S) er 


0.36 


1.61 


+0.22 



B Mixing between a CP-odd Higgs and pseudoscalar resonances 
of bottomonium 

The mixing between Higgs and resonances is described by the introduction of off-diagonal 
elements denoted by bm 2 in the mass matrix. In our case, 



VbO 

,2 



where the subindex '0' indicates unmixed states. The off-diagonal element 8m can be com- 
puted within the framework of a nonrelativistic quark potential model. For the pseudoscalar 
case under study, one can write [14] 

w = 6 (w) |fl " (0)l (B ' 2) 

Notice that bm 2 is proportional to ^, i.e. tan/5 in the 2HDM(II); high values of the 
latter implies that mixing effects can be important over a large mass region. Substituting 
numerical values (for the radial wave function at the origin we used the potential model 
estimate from [51] (^(O)] 2 = 6.5 GeV 3 ) one finds (in GeV 2 units) 

5m 2 ~ 0.146 x £ 6 (B.3) 



It is convenient to introduce the complex quantity 

'1 ] 1 

-(m 2 o - m 2 Vbo - im A f A o + im Vb0 T VbQ ) 2 + (5m 2 ) 2 



(B.4) 



19 



and the mixing quantity sin26 l = 8m 2 / A 2 , where 9 is the (complex) mixing angle of the 
unmixed rjbo resonance and A® Higgs boson giving rise to the physical eigenstates. The 
masses and decay widths of the mixed (physical) states are thus 

"4,2 - ^1,2^,2 = ^(m 2 o + m 2 Vbo - im Vb0 T Vb0 - im A oT A o ± A 2 ) (B.5) 

where subscripts 1, 2 refer to a Higgs-like state and a resonance state respectively, if m A o > 
m Vb0 ; the converse if m A o < m Vb0 . 




Figure 2: One-loop process within the SM potentially contributing to the dependence of the 
leptonic BF on the mass of the final-state leptons. 



C rjb —> £ + £ : one-loop calculation within the Standard Model 

In this appendix we consider within the SM an alternative possibility to the Higgs conjecture 
of a rising leptonic BF of bottomonium with the leptonic mass, based on the electromag- 
netic decay into a lepton pair of the rjb state subsequent to the magnetic dipole transition 
advocated in this work (see Fig. 2). Since the width for the r] b decay into two photons has 
been recently calculated elsewhere [52], obtaining the range V(j]b — > 77) — 0.4 — 0.5 keV, 
we can use the following estimate according to the SM: 

where A = (1 — 4m 2 /m 2 6 ) 1//2 . The above equation corresponds to the unitary bound due 
to the absorptive contribution of the two-photon exchange. The last inequality is readily 
obtained by setting the numerical values for tan (5 used in this work to get T Vb ^u for any 
leptonic species. 
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